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We have investigated the superfluid transition of ^He in different samples of silica aerogel. By 
comparing new measurements on a 99.5% sample with previous observations on the behaviour of 
■^He in 98% porous aerogel we have found evidence for universal behaviour of ^He in aerogel. We 
relate both the transition temperature and superfluid density to the correlation length of the aerogel. 
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The properties of bulk ■^He are well understood. The 
extreme purity of '^He at low temperatures makes it an 
ideal system to study the agreement between theoreti- 
cal and experimental results on non-conventional Cooper 
pairing in the absence of disorder. Disorder plays a cru- 
cial role in suppressing the pairing interaction in high 
Tc superconductors, the other well established non s- 
wave paired system. The superfluid transition of '^He 
confined to a sample of very porous silica aerogel was 
first reported four years ago |l|,^. The aerogel pro- 
vides a structural disorder background Q to the liquid. 
•^He is compressible, and the density can be continuously 
tuned by «30% while maintaining a fixed disorder. The 
•^He zero temperature coherence length ^O: defined as 
^0 = TiVf /ksTc, varies from 180A to over 700A as a 
function of density. Because the Cooper pairs in ^He form 
in a p-wave state, quasiparticle scattering from the aero- 
gel strands is pairbreaking ^. Thus the ^ He- in- aerogel 
system is well suited to the exploration of the effect of 
impurity scattering and disorder on the superfiuid tran- 
sition and phase diagram. 

The superfluidity of '^He in silica aerogel has been stud- 
ied using torsional oscillators ||i|,|-||, NMR (|,|-|l and 
sound propagation [ |lOyil[| techniques. These measure- 
ments show that both the superfluid transition temper- 
ature (Tc) and superfluid density (ps) of the '^He are 
suppressed by the disorder, but that the transition re- 
mains sharp pi . This suppression is sensitive to both the 
density and the microstructure of the aerogel sample. 

The simplest model for the effect of impurity scatter- 
ing on the ^He superfluid transition is the homogeneous 
scattering model (HSM) which is based on the Abrikosov- 
Gorkov model for a superconductor with magnetic im- 
purities that induce pair-breaking via spinflip scattering 
|l2| . This mechanism is similar to that of diffuse scatter- 
ing of Cooper paired ^He from a surface pl| , and is un- 
able to explain the observed behaviour. Speciflcally, the 
observed suppression of the superfluid density is much 
greater than predicted by this model. More sophisticated 
models, such as the isotropic inhomogeneous scattering 
model (IISM) |l^-[lQ| proposed by Thuneberg and co- 
workers are able to quantitatively predict the superfluid 
transition temperature of '^He in aerogel (for small sup- 
pressions) and have had success at qualitatively explain- 



ing the observed superfluid densities. 

In this Letter we present data from several different 
experiments on ^He in aerogel, including new results on 
•^He confined to a 99.5% porosity sample. This sample 
is a factor of four more dilute than any previously inves- 



tigated 1 17 1 and is crucial for understanding the evolu- 
tion from bulk '^He to a regime where impurity scatter- 
ing dominates. In comparing these different samples we 
find evidence that the relation between superfluid den- 
sity and the superfluid transition temperature of ^He in 
aerogel follows a universal behaviour, independent of the 
aerogel sample. This is significant because both of these 
quantities are individually sensitive to the microstructure 
of the aerogel, and vary greatly from sample to sample. 
We also present evidence that the suppression of Tc can 
be related to the correlation length (fa) of the aerogel 
sample. 

The aerogels used in the experiments discussed in this 
Letter were grown under basic conditions |18-p^. Un- 
der these conditions gelation is the result of diffusion 
limited aggregation of small (w30A diameter) primary 
silica particles. The aerogels are characterized by a frac- 
tal dimension (D/) related to the real space correlations, 
and a long length scale cutoff to these correlations (fa) 
above which the sample appears homogeneous. The frac- 
tal exponent depends only on the gelation process, while 
the cutoff length is also dependent on the average den- 
sity. Simulations based on the diffusion-limited cluster- 
cluster aggregation (DLCA) algorithm predict that the 
fractal exponent should lie between 1.7 and 1.9, which 
is in good agreement with small-angle X-ray scattering 
(SAXS) measurements (Table |). We note that fa in the 
most dilute sample, D, could not be inferred from the 
data as the SAXS did not extend to sufficiently small q. 





Porosity 


D/ 


ia 


Sample A 


0.98 


1.9 


1300A 


Sample B 


0.98 


1.8 


« 900A 


Sample C 


0.98 


1.8 


840A 


Sample D 


0.995 


1.7 


>2000A 



TABLE I. Parameters that characterize aerogel samples. 



Samples A and C have the same density, but were made 
under different conditions and have shghtly different frac- 
tal dimension and significantly different £,a- Both samples 
A and C have been studied with SAXS. Sample B was 
made under conditions very similar to sample C. We do 
not have direct information on its microstructure from 
SAXS, but assume here that samples B and C are essen- 
tially identical. The correlation length for sample D is 
within the range obtained from simulations based on the 
DLCA algorithm and is consistent with the SAXS data 
of Mulders on a 99.6% sample (^a >2000A). For a more 
extensive discussion the reader is referred to reference [Q . 

The 99.5% porosity aerogel used for our experiment 
was grown inside the (on average) 100 /im large pores 
of a coarse silver sinter. Previous torsional oscillator ex- 
periments have been affected by the presence of spurious 
resonances resulting from composite modes of ^He and 
aerogel whose frequency crosses the resonant frequency 
of the cell ||l|,g] . The strength of these resonances grow as 
the porosity of the aerogel sample increases, affecting the 
quality of data ||l|,||,^ . In our cell the aerogel is clamped 
to the silver sinter, and the effect of spurious resonances 
is strongly reduced. 

We operated our torsional oscillator in self-resonant 
mode near 483 Hz. The temperature in the cell was mea- 
sured using a lanthanum diluted cerous magnesium ni- 
trate ac susceptibility thermometer, thermally connected 
to the sample through a shared reservoir of ^He. Data 
was collected while the temperature increased at a rate of 
20 /iK per hour. The period shift of the oscillator as the 
superfluid '^He decoupled from the torsion head provided 
both the transition temperature of the ^He in aerogel and 
the superfluid density. 
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FIG. 1. Tlie pressure is plotted versus Tc for different 
aerogel samples. The legend is as follows: filled squares are 
cell A (98% aerogel) |I|, filled triangles are cell B (98% aero- 
gel) 1^ , filled circles are for cell C (98% aerogel) ||5| and empty 
circles are cell D (99.5% aerogel). The solid line is the super- 
fluid transition curve for bulk '^He. 

Figure shows the superfluid transition temperature 
for several different aerogel samples. Three of these sam- 
ples (A, B and C) have a nominal porosity of 98%, the 
fourth one (D) is our 99.5% aerogel. All of these mea- 
surements were done by monitoring the period shift in 



torsional oscillators filled with the aerogel and ^He. The 
difference in transition temperature between samples A, 
B and C arises from differences in the microstructure 
of the aerogels (see table | and reference Q). The rela- 
tive suppression of the transition temperature (1-Tc/Tco) 
(with Tco the transition temperature in bulk ^He) is 
larger at lower pressures (larger ^o) than at higer pres- 
sures (smaller ^g)- In view of the fact that all the aerogels 
used in these experiments have a similar fractal dimen- 
sion and primary particle size, one would expect them to 
be mutually self similar on length scales shorter than the 
fractal cutoff length. Thus for temperature dependent 
■^He coherence lengths £,{T) shorter than ^q, the ratio of 
density of silica sampled at two different £,{T) is inde- 
pendent of aerogel density. Experimentally one observes 
a strong dependence of the suppression of T^ on the aero- 
gel density. 
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FIG. 2. The superfluid transition temperature for '^He in 
aerogel, scaled by the bulk transition temperature, is plotted 
for several different samples as a function of ^o . The symbols 
are the same as shown in Fig. |l[ 

In figure we show the dependence of the relative sup- 
pression of Tc as a function of ^o- At short coherence 
lengths, the relative suppression is small for all sam- 
ples. However, when the pressure is reduced, increasing 
^0, this suppression shows a marked dependence on the 
microstructure of the aerogel. The transition tempera- 
ture shows an evolution from a strong impurity scatter- 
ing regime when the '^He is confined to a 98% porosity 
aerogel sample towards the behaviour of bulk '^He in the 
99.5% sample. 

Motivated by figure g we plot the relative suppression 
of Tc against ^o scaled by the aerogel correlation length 
^Q. The aerogel-limited mean free path (Ig) would be an- 
other natural choice to compare against ^o- However, this 
length scale has not been independently measured. If we 
use the HSM to determine Ig from Tc/Tco, the values of Ig 
show a strong pressure dependence 122] . Figure shows 
that the scaled transition temperature depends solely on 
the ratio of ^o to ^^ independent of the aerogel density. 
The error bars for sample D result from the high and low 
estimates (2000A and 3000A respectively) for the cor- 
relation length for a 99.5% porosity sample from DLCA 
simulations. 
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FIG. 3. Relative suppression of Tc versus ^o/Ca- The cor- 
relation length for the 99.5% sample is taken to be 2200A with 
the error bars corresponding to the high and low values of 
3000 A and 2000 A. 

We determined the temperature-dependent superfluid 

density of ^He in 99.5% aerogel using the shift in res- 
onant frequency of our torsional oscillator upon warm- 
ing. A small, temperature dependent contribution due 
to bulk ^He in the cell was subtracted. The remaining 
shift was scaled by the period change due to filling the cell 
with '^He at 50 mK and by the tortuosity (measured with 
''He) to obtain the superfluid density. Figure |4| shows the 
bare superfluid density (p^/p) plotted versus the reduced 
temperature T/T^ for different aerogel samples and pres- 
sures. The bare superfluid density is obtained from ps by 
stripping away the Fermi liquid factor [p3| according to: 



1 i^+kFin 



1 



3 ^ p 



(1) 



and is equivalent to l-Y(T) where Y(T) is the temper- 
ature dependent Yosida function for bulk ^He. 




FIG. 4. We plot the bare superfluid density versus re- 
duced temperature for cell A (98%), flUed symbols, and cell 
D (99.5%), empty symbols. 

As with the transition temperature T,,, the Pg for these 
two different samples are similar at high pressures and 
are both close to the bulk value. At the lowest pressure, 
there is a factor of five difference in p^/p at the same 
reduced temperature between the 99.5% sample and the 
98% sample; p^/p is more strongly suppressed by the 



aerogel than Tc. There is also a large difference between 
the suppression factors of Tc and p^ in the same aero- 
gel sample. This large suppression in ps with respect 
to Tc is consistent with measurements made on other 
samples [EH| , and cannot be explained with homogenous 
scattering models for '^He in aerogel [Q . 
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FIG. 5. The bare superfluid density at 0.8Tc scaled by 
the bare bulk superfluid density is plotted versus the square 
of the scaled Tc. The symbols are the same as in Fig. hi The 
dashed line is the theoretical curve from the HSM. 

In order to better understand the superfluid transi- 
tion of '^He in aerogel, in figure we plot p'^^/ps at 
0.8Tc against (Tc/Tco)^ for several aerogel samples. The 
error bars shown for cell B arise from a spurious reso- 
nance that introduces uncertainty in the determination 
of Ps- As in figure ||, the data collapses on a universal 
curve. The dashed line is the prediction for a homoge- 
neously scattering model based on the Abrikosov-Gorkov 
equation flj] . This plot compares aerogels with different 
densities — there is a factor of four difference in the aver- 
age impurity density between the 98% samples and the 
99.5% aerogel sample. Furthermore, the coherence length 
of the Cooper pairs varies from 180 A to 600 A over this 
data, yet the very strong pressure dependence shown in 
figures pi and has been factored out in this plot. 

The correlations in the aerogel will affect the suppres- 
sion of Tc and the evolution of p^ relative to a homogeous 
disorder, but this suppression apparently depends only 
on the correlation length of the sample and possibly the 
fractal exponent. The steep slope of the data in figure H 
is evidence that the fractal nature of the aerogel plays an 
important role in the development of ps pq | . Since all of 
these samples were base-catalyzed, the behaviour of '^He 
in each aerogel is determined mainly by ^a ■ As long as ^o 
is much less than ^^ , the disorder sampled by the ensem- 
ble of Cooper pairs should be insensitive to changes in 
the temperature-dependent coherence length ^(T), until 
^(T)w ^Q. That is, the system has a conformal symmetry 
normally absent in disordered systems. The evidence of 
this one-parameter scaling is displayed in figure q p\ is 
a function of Tc/Tco (or equivalently Co/Ca) only. This 
behaviour is reminiscent of the compilation of data from 
disordered high Tc materials by Franz et al Efl] with the 
exception that the relative suppression of Tc/Tco and 



Ps/Pso in the high Tc materials do not follow a univer- 
sal behaviour, presumably because the impurities are not 
fractally correlated. 

In order to compare the behaviour of Tc and p^ of '^He 
in aerogel with "^He in bulk it will be necessary to under- 
stand precisely how the fractal disorder affects the su- 
perfluid pairing mechanism. One test for models of non- 
conventional Cooper pairing in the presence of disorder 
would be to predict the functional form of the univer- 
sal curve for p^ versus (Tc)^ for ■^He in (base-catalyzed) 
aerogel. The IISM of Thuneberg et al. predicts a re- 
lation between p\/ p\n and Tc/Tco similar to the trend 
illustrated in figure |, showing behaviour very different 
than Abrikosov-Gorkov model n,15|. This model does 
not explicitly consider the fractal nature of the aerogel, 
but shows how inhomogeneities in the disorder can lead 
to a large suppression of pi relative to T^. 

In this letter we have presented data from our mea- 
surements on ^He in a very dilute 99.5% porous silica 
aerogel. The values of Tc and p\ fall between those of 
■^He in bulk and '^He in denser aerogel samples. We also 
present Tc and p\ data for '^He in aerogel experiments 
performed at Cornell that show universal behaviour that 
can be traced to the fractal structure of the aerogel. In 
order to more fully understand this exciting physical sys- 
tem, more attention must be devoted to understanding 
the microstructure of the aerogel. Specifically, the uni- 
versal scaling discussed above depends strongly on the 
fact that the fractal exponent real-space correlations is 
similar for all the aerogel samples. It would be inter- 
esting to study neutrally catalyzed silica aerogels, which 
have a different fractal exponent than the base-catalyzed 
samples 1l9J. '^He would be expected to follow different 
universal behaviour when confined to base-catalyzed and 
neutrally catalyzed aerogels. As yet, there is no clear 
theoretical picture for why the transition temperature 
TcShould depend on the ratio ^o/Ca- The explanation for 
this behaviour will provide insight into how correlated 
disorder affects non-conventional Cooper pairing. 
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